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1.  Introduction 


The  computation  of  electric  and  magnetic  fields  at  dielectric  boimdaries 
must  take  into  account  interrelated  boimdary  conditions,  two  involving 
normal  components  of  the  electric  and  magnetic  fields  and  four  involving 
tangential  components  of  these  fields.  In  a  few  simple  cases,  it  is  possible 
to  determine  the  configuration  of  the  polarization  or  propagation  modes 
from  general  principles  and  to  apply  boundary  conditions  separately  to 
each  mode.  For  example,  in  the  derivation  of  the  Fresnel  equations,  which 
describe  refraction  of  a  plane  wave  at  a  planar  surface,  boundary  condi¬ 
tions  are  applied  separately  to  the  s  and  p  polarizations,  with  the  electric  or 
the  magnetic  fields  in  the  plane  of  the  surface.  There  are  only  three  bound- 
aty  conditions  for  each  polarization,  and  the  magnitudes  and  propagation 
directions  for  each  are  determined  readily.  Similarly,  the  lowest-order 
propagation  modes  of  a  dielectric  fiber,  which  have  no  azimuthal  field 
variations,  are  transverse  electric  and  transverse  magnetic  [1],  each  with 
only  three  components;  the  magnitudes  of  the  fields  and  propagation  con¬ 
stants  can  be  computed  separately  for  each  mode. 

When  the  modal  configurations  are  not  known,  two  approaches  can  be 
taken  to  determine  the  fields.  In  one  approach,  four  independent  field  con¬ 
figurations  that  account  for  all  possibilities  permitted  by  Maxwell's  equa¬ 
tions  are  constructed  in  such  a  way  that  the  corresponding  field  compo¬ 
nents — tangential  and  normal  electric  and  magnetic  components — on  each 
side  of  the  interface  are  proportional  everywhere  on  the  interface.  The  con¬ 
tinuity  conditions  at  the  boundary  are  applied  to  four  of  the  components 
(e.g.,  the  two  tangential  components  of  the  electric  and  magnetic  fields)  to 
determine  relative  magnitudes  of  the  four  field  configurations.  This 
approach  was  used  by  Mie  [2]  to  describe  scattering  from  a  dielectric 
sphere,  and  by  Hondros  [3]  and  Hondros  and  Debye  [4]  to  describe  propa¬ 
gation  in  a  dielectric  fiber. 

The  other  approach  uses  electromagnetic  scalar  potentials,  which  are  pairs 
of  scalar  solutions  to  the  wave  equation  from  which  the  electric  and  mag¬ 
netic  fields  can  be  derived.  Debye  [5]  used  a  pair  of  potentials  (known 
today  as  the  Debye  potentials)  to  describe  scattering  from  a  dielectric 
sphere,  which  is  the  same  problem  that  was  solved  by  Mie.  Debye  devel¬ 
oped  expressions  for  the  fields  in  terms  of  these  functions,  and  he 
described  boundary  conditions  that  the  scalar  functions  and  their  first 
derivatives  were  to  satisfy  on  the  interface  (see  eq  (68)).  To  satisfy  these 
conditions,  Debye  found  scalar  functions  in  which  the  values  of  the  two 
functions  and  their  normal  derivatives  were  proportional  on  each  side  of 
the  interface. 

In  addition  to  the  Debye  potentials,  other  pairs  of  electromagnetic  scalar 
potentials  have  been  investigated.  Whittaker  [6],  who  appears  to  have 
been  the  first  person  to  realize  that  electromagnetic  fields  could  be  derived 
from  two  scalar  potentials,  described  a  pair  from  which  fields  due  to  mov¬ 
ing  point  charges  could  be  computed.  Green  and  Wolf  [7]  used  a  different 
pair  of  potentials  to  compute  the  energy  and  momentum  of  electromag¬ 
netic  fields  and  to  construct  a  Lagrangian.  Nisbet  [8]  developed  a  general 
transformation  that  relates  possible  pairs  of  scalar  potentials,  and  Stratton 
[9]  gave  a  simple  expression  for  magnetic  vector  potentials  in  terms  of  an 
arbitrary  set  of  scalar  potentials  (see  eq  (66)). 
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More  recently,  a  representation  of  electromagnetic  fields  in  terms  of  a  dif¬ 
ferent  pair  of  scalar  functions  was  developed  by  Morgan,  et  al  [10]  and 
extended  by  Morgan  [11]  for  computing  propagation  through  nonuniform 
dielectrics.  These  scalar  functions — coupled  azimuthal  potentials — ^are  not 
electromagnetic  scalar  potentials  in  the  sense  discussed  above,  because 
they  do  not  independently  satisfy  the  wave  equation.  Instead,  they  are 
solutions  to  a  pair  of  coupled,  second-order  differential  equations.  Solu¬ 
tions  are  found  by  assuming  Dirichlet  boxmdary  conditions  at  a  geometric 
boimdary  of  a  body  and  computed  across  a  dielectric  interface  somewhere 
within  the  body  by  a  finite-difference  method  that  begins  at  the  geometric 
boimdary.  Moser,  et  al  [12]  used  coupled  azimuthal  potentials  to  calculate 
frequencies  of  scattering  resonances  of  bodies  of  revolution  for  application 
to  target  recognition.  Several  other  applications  are  referenced  in  the  paper 
by  Morgan  [11]. 

The  requirement  that  both  the  scalar  functions  and  their  derivatives  be 
proportional  everywhere  on  each  side  of  a  dielectric  interface,  such  as  the 
arrangement  represented  for  the  Debye  potentials  in  equation  (68),  is  diffi¬ 
cult  to  satisfy  for  complex  geometries.  Because  of  the  spherical  symmetry 
of  his  scattering  object,  Debye  was  able  to  express  his  scalar  potentials  in 
terms  of  spherical  wave  functions,  and  to  satisfy  the  four  continuity  condi¬ 
tions  with  those  functions.  Most  often,  however,  the  forms  of  the  functions 
are  not  known  from  general  considerations,  and  a  search  for  functions  that 
match  at  the  interface  can  be  difficult. 

In  this  report,  it  is  shown  that  the  scalar  function  that  is  zero  at  the  surface 
and  the  one  whose  normal  derivative  is  zero  can  be  used  as  scalar  poten¬ 
tials.  Using  these  functions  reduces  the  number  of  proportionality  condi¬ 
tions  at  the  boundary  that  must  be  considered  from  four  to  two,  which 
relates  the  normal  derivative  of  the  first  function  and  the  value  of  the  sec¬ 
ond  function  on  each  side  of  the  interface.  Boundary  conditions  at  the 
interface  are  expressed  in  terms  of  these  fimctions,  and  expressions  for  the 
electric  and  magnetic  fields  in  terms  of  these  functions  are  so  derived. 

It  should  be  pointed  out  that  the  scalar  potentials  that  satisfy  the  Dirichlet 
and  Neumann  conditions  are  used  implicitly  in  the  solution  of  boundary 
value  problems  involving  perfect  conductors.  Here,  a  tangential  compo¬ 
nent  of  the  electric  field  is  equated  with  a  scalar  function  that  is  zero  on  the 
conducting  surface,  and  a  similar  tangential  component  of  the  magnetic 
field  is  equated  with  a  scalar  function  whose  normal  derivative  is  zero.  For 
example,  in  the  determination  of  fields  in  a  conducting  waveguide,  the 
longitudinal  electric  and  magnetic  fields  can  be  represented  by  the  two 
scalar  functions,  and  the  transverse  fields  can  be  derived  from  these 
functions. 

In  the  following  sections,  a  set  of  equations  is  derived  that  describes 
boundary  conditions  in  terms  of  magnetic  vector  potentials.  The  vector 
potentials  are  expressed  in  terms  of  scalar  functions  that  obey  the  Dirichlet 
and  Neumann  boundary  conditions.  Solutions  in  terms  of  the  scalar  func¬ 
tions  are  found  for  the  equations  that  describe  boundary  conditions  at  a 
dielectric  interface.  These  solutions  are  combined  in  a  way  that  satisfies  the 


radiation  condition.  There  are  two  possible  combinations,  which  are  iden¬ 
tified  with  propagation  modes  or  polarizations.  Through  the  equations 
that  relate  electric  and  magnetic  fields  to  the  vector  potentials,  these  fields 
are  also  expressed  in  terms  of  the  scalar  functions.  To  demonstrate  this 
approach,  the  methodology  is  used  to  derive  the  Fresnel  equations.  The 
methodology  is  also  applied  to  find  fields  bounded  by  a  conducting 
surface  where,  unlike  the  dielectric  interface,  there  is  no  continuity 
requirement  at  the  boimdary,  and  differences  between  the  two  situations 
are  discussed.  Finally,  the  scalar  potentials  defined  in  this  paper  are  com¬ 
pared  with  those  used  by  Debye  [5]  and  by  Morgan,  et  al  [10],  and  similari¬ 
ties  and  differences  are  discussed. 

A  two-dimensional  approach  is  taken  here  by  representing  fields  and 
potentials  as  products  of  functions  of  two  coordinates  and  the  exponential 
exp(-ifijjz),  where  z  is  the  third  coordinate.  The  results  are  directly  appli¬ 
cable  to  a  uniform  waveguide  (where  the  third  coordinate  is  the  longitudi¬ 
nal  distance),  or  to  a  rotationally  symmetric  scatterer  (where  the  third 
coordinate  is  in  the  azimuthal  direction  about  the  axis  of  symmetry). 
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2.  Boundary  Conditions  for  Vector  Potentials 


The  boundary  conditions  for  the  vector  potentials  arise  in  a  natural  way  by 
multiplying  them  by  step  functions  that  change  from  zero  to  one  at  the 
interface,  and  taking  derivatives  in  the  usual  way.  The  coefficients  of  the 
delta  function  that  result  from  the  differentiation  are  the  boundary  condi¬ 
tions,  and  the  coefficients  of  the  step  function  make  up  the  homogeneous 
wave  equation  for  regions  of  constant  permittivity  and  permeability.  This 
approach  was  demonstrated  for  the  electric  and  magnetic  fields  by  von 
Roos  [13]  and  is  applied  here,  in  a  similar  way,  to  sets  of  vector  potentials 
that  terminate  at  a  conducting  surface,  and  which  are  transmitted  through 
a  dielectric  interface. 

Consider  first  the  fields  produced  on  a  perfectly  conducting  surface. 
Assume  initially  that  the  surface  is  planar,  with  a  coordinate  system 
shown  in  figure  1.  The  normal  to  the  surface  is  in  the  x-direction.  In  terms 
of  vector  potentials,  the  Maxwell's  equations  for  the  fields  are 

S7xVxA  +  eaH^  +  Solii^^O  =  ^jSix) ,  (1) 

and 

v2<i,  +  v.^  =  -lpj(x),  (2) 

where  S{x)  is  a  delta  function,  and  (A,<I>)  are  products  of  these  solutions 
and  the  step  function; 

A{x,y,z)=T(x)Ao(;x,}/,z) , 

<t>  (x,  y,  z)  =  T  (x)  Op  (x,  y,  z)  , 

where  (Ao,^o)  are  solutions  to  the  homogenous  equations  that  result  when 
the  right  sides  of  equations  (1)  and  (2)  are  zero. 

The  step  function  is  defined  as 

r(x)  =  0  ifx<0  , 

=  1  ifx>0.  « 


Figure  1.  Directions  of 
radiation  incident  on 
a  conducting  surface. 


0  jc 


8 


Substituting  equation  (3)  into  equation  (1)  and  applying  vector  identities 
gives 


•Ao  +  e„/i 


t{x) 


+ 


+ 


’dx^ 

d^T 

'dx^ 


+  2 

+  2 


^yo  dr 
dx  dx 


^xodt 

dy  dx 


y 


dr  ^xo  dr 
Bx  dx  dz  dx 


-MjS{x) , 


(5) 


where  Ayg  is  the  y-component  of  Ag  and,  similarly,  for  the  other  coordi¬ 
nates.  In  like  manner,  equation  (2)  may  be  expressed  as 


/  2  B 

ao; 

1  “  se  it 

V-Ag+EgH 

r{x) 


O  — +  2^^— + 
°  dx^ 


1 


=  -  pS{x)  . 


(6) 


The  first  lines  of  equations  (5)  and  (6)  are  identically  zero,  because  the  com¬ 
ponents  of  (AcOp)  satisfy  the  wave  equation  and,  together,  they  obey  the 
Lorentz  condition 


BAxo  ^■'^yo  1  ^^0 


The  other  terms  of  equations  (5)  and  (6)  define  the  boundary  conditions. 
To  evaluate  them,  use  the  derivatives  of  t,  given  by  [14]  as 


These  terms  become 


dh^  S{x) 
dx^~  X  ‘ 


S(^) 


A 

X 


w  _  2  ^  ^^xo 


Bx  By 


y  +  5{x) 


A 

X 


TO  2  ^20  ^^XO 


Bx  Bz 


z  = 


^ix)^^{iyy+izz) 

and 


(8) 


(9) 


3tc 


BAxo 

Bt 


(10) 
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Terms  with  x  in  the  denominator  have  meaning  only  if  the  numerators  are 
zero  at  the  boundary.  In  this  case,  the  ratio  can  be  expressed  by  a  deriva¬ 
tive.  With  these  substitutions,  the  boimdary  conditions  become 


and 


<^0  =  0, 
^ZO  “  0  , 

^yo  “  0  / 


^xo 

3tr 

dt 

dA„o 

dx 

az 

dAyo 

^xo 

dx 

dy 

=  -Hjy  . 


(11) 


(12) 


The  corresponding  boundary  conditions  at  a  dielectric  interface  are 
derived  by  setting  the  right  sides  of  equations  (1)  and  (2)  to  zero  and  repre¬ 
senting  the  potentials  and  permittivity  by 


0  =  +  r<t>2  , 

A  =  {\-i)A-^  +  zA2  , 
e=(l-'^ei  +  Te2  , 

where  (Ai,<I>i)  and  (^2,02)  are  the  solutions  to  the  homogeneous  wave 
equations  in  media  1  and  2,  and  Ei  and  £2  are  the  permittivities  in  the  two 
media.  Similar  to  the  potentials  described  above,  (Ai,Oi)  and  (A2,02)  obey 
the  Lorentz  condition.  Applying  vector  identities  as  before  and  gathering 
coefficients  of  the  delta  functions  gives  the  boundary  conditions  at  the 
dielectric  interface  in  terms  of  vector  potentials 


0i=02  , 

^y\  -  ^yl  > 


aA^i 

302 

dAx2 

dx 

dt 

-^2 

dx 

dt 

dA^i 

dA^ 

dA^ 

dx 

ofe 

dx 

dz 

aAyl 

aA^i 

dAy2 

dA^ 

dx 

dy  \ 

[  ^ 

(14) 


(15) 


These  boundary  conditions  for  the  vector  potentials  are  equivalent  to  the 
usual  boundary  conditions  for  electric  and  magnetic  fields.  Equation  (14)  is 
related  to  the  continuity  of  the  tangential  components  of  the  electric  field  E 
and  the  normal  component  of  the  magnetic  induction  B  across  the  inter¬ 
face.  Members  of  equation  (15)  are  statements  of  the  continuity  of  the  nor¬ 
mal  component  of  the  electric  displacement  D  and  the  tangential  compo¬ 
nents  of  the  magnetic  intensity  H  at  the  interface. 


3.  Scalar  Functions  at  a  Boundary  Surface 

Two  sets  of  potentials  will  be  computed  separately.  In  the  first  set,  three  of 
the  potentials  are  zero  at  the  boundary  and  the  normal  derivative  of  the 
vector  potential  component  normal  to  the  surface  is  zero.  These  are  the 
conditions  of  equation  (11),  and  the  potentials  in  this  group  are  also  the 
expressions  for  the  potentials  at  a  perfectly  conducting  interface.  The  sec¬ 
ond  group  of  solutions  obeys  a  set  of  conjugate  boxmdary  conditions  for 
which  the  normal  derivatives  of  the  first  three  potentials  and  the  fourth 
potential  itself  is  zero.  We  express  each  of  the  vector  potentials  in  terms  of 
the  functions 


M(x,y)exp[-i(©2Z-®f)]  , 
n(x,y)exp[-i(a)2  2-tof)]  , 

where  u  (x,y)  and  the  normal  derivative  of  v  (x,y)  are  zero  on  the  boundary. 
Anticipating  that  we  will  want  expressions  for  fields  that  are  valid  for  sur¬ 
faces  of  any  orientation  of  the  boimdary  in  the  (x,y)  plane,  we  seek  repre¬ 
sentations  for  vector  potentials  that  are  rotationally  invariant  in  this  plane. 
A  set  of  expressions  for  the  first  group  of  solutions  that  satisfies  these  con¬ 
ditions  is: 


^zl  > 

du-i  dvj^ 

_  auj  .  doi 


(17) 


where  (0ii,  a^ii,  dn,  an)  are  constants.  In  these  equations,  the  subscripts  on 
the  potentials  (Oi,Ai)  and  on  the  scalar  functions  {ui,Vi)  refer  to  region  1. 
The  constants  have  two  numerical  subscripts:  the  first  subscript  refers  to 
the  region,  and  the  second  subscript  to  the  first  or  second  set  of  solutions. 


The  coefficients  (^n,  0^1,  dn,  flu)  are  evaluated  by  substituting  equation 
(17)  in  which  they  are  regarded  as  variables  into  equation  (15)  to  produce 
sets  of  linear  equations.  The  Lorentz  condition  is  used  to  eliminate  d^.  In 
matrix  form,  the  equations  for  the  remaining  coefficients  are 


K^Col 

(0^  0) 

0)Kiq^ 

hi 

li 

kf 

dui 

RTl  (O2  (O/C^ 

^zll 

dx 

kl 

0 

0 

k\v^ 

hu-^/dx 

«11 

1^3 

(18) 


11 


for  region  1  and,  similarly,  for  region  2.  Here  k;  and  are  the  specific  per¬ 
mittivity  and  the  square  of  the  magnitude  of  the  propagation  vector,  given 
by 


K-,  =  e,/£o  ' 

r  2  2 

ki  =  Krj-o)^  , 
c 


(19) 


'll  “  ^ — 7^  •  assume  qi  is  constant  on  both  sides  of  the  boundary. 

oti'^f  doc 

Without  loss  of  generality,  we  can  normalize  {u,v)  so  that  their  derivatives 
are  equal  at  the  boundary 


Buj  du2 
dx  dx 

do  2  (20) 

As  a  consequence  of  this  normalization,  v-y  =  V2  and  q^  =  q2  hold. 

The  only  quantities  in  equation  (18)  that  are  characteristic  of  the  medium 
are  k:,  and  ki^.  Equation  (18)  can  be  satisfied  by  solutions  that  lead  to 
(^1/^2/ ^3)  that  are  not  functions  of  either  of  these  quantities.  A  set  of  solu¬ 
tions  for  which  this  is  true  is 


(O/C 

{l/K,)Q}^q 

^Zll 

= 

(ojc 

/ 

(Olc^ 

f 

0 

0 

0 

{oxo,)/{ck,]^ 

(21) 


where  i  is  1  or  2.  The  first  of  these  solutions  defines  gauge  transformations 
and  the  other  two  lead  to  nonzero  fields. 


The  second  group  of  solutions  to  equation  (15)  can  be  written  in  terms  of 
(u,v)  as 

Oi  =  012^1  ^ 

3^1  dui 

do^  du-i 


(22) 


for  region  1  and,  similarly,  for  region  2.  As  in  equation  (17),  these  solutions 
are  rotationally  invariant  in  the  {x,y)  plane.  Their  substitution  into  equa¬ 
tion  (15)  leads  to  the  following  equations  for  the  constants 


dui  9^2 


=  £2 


dv 


du 


1 


^Zl2~  - (0yd 


dx 


'2^12' 


dx^ 


+  d 


dvi 

d^v 


9^2  3^2  ^2 

<l,^—-(oa22-^-cod^— 


422' 


dV2 

dx 


du 


(OyU 


1 


2  ,  ^2 

-  a,d^ 


=  a 


22' 


9r^ 


+  d 


d^V-y 


■21 


dxdy 


(23) 


Both  sides  of  each  of  these  equations  are  identically  zero  because  the 
derivatives  of  u  and  v  that  appear  in  these  equations  are  zero  at  the  bound¬ 
ary.  The  boundary  conditions  represented  in  equation  (14),  which  were 
identically  zero  for  the  potentials  of  equation  (17),  are  used  to  evaluate  the 
constants  in  equation  (22).  Substitution  of  equation  (22)  into  the  first  three 
members  of  equation  (14)  gives: 


^12^1  -  <(>22^2  ' 

®  212^1  =<^222^2  ' 

9^1  9ui  9^2  du2  (24) 

The  Lorentz  condition  may  be  used  to  eliminate  one  of  the  constants,  as 
before.  The  resulting  equations  can  be  expressed  in  matrix  form  as 


1 

0 

0 

^i2 

w, 

0 

1 

0 

^zi2 

W2 

®2^ 

1 

^i2 

^3 

kf 

As  with  equation  (18),  equation  (25)  may  be  satisfied  by  sets  of  constants 
that  yield  values  of  (Wi,W2,W3)  that  are  not  fimctions  of  k;  or  k,^.  A  set  of 
linearly  independent  solutions  is 


^i2 

a)/c 

0 

0 

^*212 

= 

(ojc 

/ 

(Olc^ 

/ 

0 

0 

(<y/c^) 

L^l 

Mi2 

. 

«2 

(26) 
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As  with  equation  (21),  the  first  of  these  solutions  defines  gauge  transfor¬ 
mations  and  the  others  define  propagation  modes. 


The  electric  and  magnetic  fields  can  be  computed  by  substituting  the 
potentials  described  by  equations  (17)  and  (22)  into  the  equations  for  the 
electric  and  magnetic  fields,  which  becomes 


B=VxA , 


(27) 


It  is  convenient  to  define  the  quantities 

=  (l Af)  (^,1  -  Oflzfl)  Ui  +  {(0,  (t>i2  -  v]  , 

where  the  subscript  i  denotes  regions  1  or  2.  In  terms  of  these  quantities, 
the  electric  and  magnetic  fields  are 


(28) 


'XI 


dbi 

K^eodei 

ab,  dei 

t-  ,  —  />■> _  4.  _ 

'dx 

^  ' 

"dr  ' 

dbi 

_  4. 

KjCOdCi 

db;  dc; 

r  _  _  /•,, '  j.  : 

% 

,2  dx  ' 

1 

1 

■bi, 

E^i  =  ikfei  . 

(29) 


Substitution  of  the  first  solutions  of  equations  (21)  and  (26)  into  equation 
(29)  yields  values  of  fe,  and  Ci  that  are  zero,  which  verifies  that  these  solu¬ 
tions  define  gauge  transformations. 

General  expressions  for  the  sets  of  constants  in  equations  (21)  and  (26)  that 
lead  to  nonzero  fields  may  be  constructed  as  linear  combinations  of  the 
second  and  third  solutions: 


and 


(1//C,)  ftj. 

(l//q)  (0^ 

^izl 

-^1 

(Olc^ 

0 

0 

{a)co,)/{ck^f 

0 

0 

^iz2 

=  h 

(O/c 

+  f4 

0 

‘^il 

{(O/c^) 

fl+— 1 

co/c^ 

0), 

(Oz 

(30) 


(31) 


where  (ti,  t2,  ^3,  t^)  are  constants.  Substitution  of  equations  (30)  and  (31) 
into  equation  (28)  leads  to  expressions  for  fc,  and  e,  in  terms  of  (fj,  t2,  ti): 


bi  =  - 


(O  1 

P-  0), 


1  + 


V 


+  t/ 


^,  =  - 


K 


h-h4 

k} 


'■i  y 


Wi¬ 


fi)  ®2 


1  (0^ 
kfp 


h^i  ■ 


(32) 


Because  the  values  of  (b„  e,)  are  independent  of  the  orientation  of  the  coor¬ 
dinate  system  in  which  (u„  i;,)  were  defined,  equation  (32)  holds  for  any 
point  on  the  interface,  provided  that  q  has  the  same  value  at  every  point; 
the  functions  u  and  v  must  be  divided  into  components  for  which  this  is 
the  case.  For  example,  if  the  interface  is  planar,  u  and  v  might  be 


The  ratio  q  would  be 


M~w(x)exp(-ifi)yi/)  , 
o  ~  u  (x)  exp  (- i  fijy  y)  . 


(33) 


-ia)yV(x) 

aM(x)/atc  ■ 

The  ratio  y  is  a  function  of  (Oy,  implying  that  propagation  across  the  inter¬ 
face  would  be  different  for  components  of  different  (Oy.  The  condition  that 
q  be  constant  and  the  conditions  of  equation  (20)  are  requirements  for  the 
functions  u  and  v  at  the  interface. 

It  was  assumed  initially  that  the  boundary  of  the  region  where  fields  are 
calculated  was  planar.  Since  the  fields  in  this  two-dimensional  problem  are 
proportional  to  u  and  v  or  their  first  derivatives,  this  restriction  can  be 
dropped  because  surface  curvature  enters  only  in  the  computation  of  sec¬ 
ond  derivatives. 
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4.  Radiation  Condition 


In  addition  to  the  boundary  conditions  at  the  dielectric  interface,  another 
condition  must  be  satisfied:  radiation  that  propagates  across  the  interface 
must  be  outgoing  or  exponentially  deca)ing.  Iri  the  configuration  of  figure 
2,  radiation  is  incident  from  region  1  and  radiation  outgoing  in  region  2. 
This  configuration  places  constraints  on  possible  values  of  (fj,  t2,  U).  To 

calculate  these  constants,  we  can  express  U2  and  V2  as  functions  of  outgoing 
and  incoming  radiation  functions /and  g  as 

!’2  +  A/+fe. 

where /represents  outgoing  or  exponentially  decaying  radiation  and^ 
represents  incoming  or  exponentially  growing  radiation.  The  quantities  i»2 
and  02  become 


and 


7  LU  ± 

b-,=-{a„— — 


kl 


kl  r ' 


(36) 


-f- 


(olq 


^4 


C 


^  4 


-B— 

kl 


1  (0^ 


h}g 


(37) 


Figure  2.  Local 
coordinate  system  on 
interface  between  two 
dielectrics. 


region  2 
t  =  t2 
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The  coefficients  of  ^  must  be  zero.  Equations  for  a  and  P  are 
0)  1  \  (  1  CO 


1  0)h  1  0)^ 


Because  these  equations  are  homogeneous,  they  can  only  have  nonzero 
solutions  for  a  and  P  if  the  coefficients  of  oc  and  P  in  each  equation  are  pro¬ 
portional.  Let  T  be  the  constant  of  proportionality,  defined  by 


c 


In  terms  of  T,  equation  (38)  is 


Of  ^3  1  H - —  +  ^4  4-  pr 


,  KsltoVc^) 


^3  —  0  , 


or  - ^ - Tts  +p - - - 


The  coefficients  of  a  must  also  be  proportional,  as  in 
(  ^z\  K2{o)^/c^)q 

I  - 

Equations  (39)  and  (41)  may  be  used  to  eliminate  f2  and  from  the  expres¬ 
sions  for  hi  and  ei.  Equation  (40)  may  be  used  to  express  fj  in  terms  of  the 
ratio  p/ a.  These  expressions  for  (ti,  ti,  h,  t^)  may  be  substituted  into 
expressions  for  and  Ci,  which  are  similar  to  equations  (36)  and  (37),  with 
ki  replaced  by  k^.  The  result  is 

,  (o  1  \J  Ka]  ^f^2p\  „  dco/c^] 

A  similar  expression  for  medium  2  is 

,  ^0)  1  fCifp  \ 

&2  =  T— - 5-  -  “2-^2  / 

kl  l«  } 

„  1  ^  (43) 

^2  =  -^  -“2-^2  • 


n 


All  the  electric  and  magnetic  fields  may  be  derived  from  {b,  e)  using  equa¬ 
tion  (29).  The  parameter  T  defines  propagation  modes.  It  is  often  conve¬ 
nient  to  define  modes  with  T  =  0  and  T  =  <»,  but  it  is  not  necessary  to  do 
this.  Any  two  values  of  T  will  define  linearly  independent  sets  of  (b,  e)  in 
terms  that  may  express  all  other  possibilities. 

Although  this  derivation  was  carried  out  with  a  radiation  condition  in 
mind,  the  result  is  more  general.  The  procedure  is  valid  because  the  two 
independent  solutions  (u,  v)  to  the  wave  equation  can  be  expressed  as  a 
linear  combination  of  two  different  solutions  (f,  g).  The  functions  (f,  g)  may 
represent  solutions  that  are  physically  admissible  or  inadmissible.  In  this 
case,  the  admissible  solutions  were  outgoing  or  exponentially  decaying, 
and  the  inadmissible  solutions  were  incoming  or  exponentially  growing. 


5.  Fresnel  Equations 


The  simplest  interface  between  two  dielectric  media  is  a  plane.  Figure  3 
shows  incident  and  reflected  waves  in  medium  1,  and  an  outgoing 
refracted  wave  in  medium  2.  We  assume  that  the  waves  are  planar,  with 
directions  of  propagation  in  the  x-z  plane.  The  radiation  has  no  depen¬ 
dence  on  the  y-coordinate.  In  medium  2,  the  functions  u  and  v  are 


«2  =  X  exp  (-  sin  x  , 

^2  =  exp  (- cos  X  . 

For  outgoing  and  incoming  radiation,  the  functions/ and  are 

f=exp[-i[o)^z  +  (o^x)]  , 
g  =  exp[-i[(a^z-(o^x]]  . 

The  constants  a  and  p  are 

a=-i{af2)  , 

P=l/2  , 
p/a=i/a  . 

Functions  Ui  and  in  medium  1  are  defined  so  that 

dui/^x  =  ^U2/dx  , 

Vi  =  V2  . 

With  these  normalizations,  Wi  and  are 

«!  =  fl  (<0x2/ ^xi]  ®xp  (-  icD^zj  sin  cd^^x  , 

=  exp  (-  ioj^z'j  cos  co^ix  . 


(44) 


(45) 


(46) 


(47) 


(48) 


Figure  3.  Planar 
interface  between  two 
dielectrics  and 
directions  of 
propagation  of 
radiation. 
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The  ratio  of  derivatives  q  is  zero,  because  v  has  no  dependence  on  y.  With 
this  value  of  q,  equatioris  (42)  and  (43)  become 


c^O)^\akl  k\  I 
i  'I  1 

e^= - -Ml  — -vi  , 

^2  fl/c?  ,  ki 


b-,  =  T- 


0)^2  l(i 


-U2-V2  , 


^2  =  71  -«2-^2  . 


Substitution  of  equations  (44)  and  (48)  into  these  expressions  gives 


7  rri  ^  ^  • 

bi  =  T——  i 


CO  >^2  .(»x2  1  .  1  1  ,  ■  \ 

ITT  ^  — TiSinoj^ia:-  — cosfo^ix  expl-My^z)  , 


\  ^xi  k^2 
.  f^2  (0x2  1  . 


ei=  t - r 

1  >^1  kl 


(0x2  1  .  1  \  f  ■  \ 

—  —  sm  <0^1  X  -  —  cos  (O^jX  exp  (-  tco^z)  , 


7  rj-,  ^  ^2  1  .  \  (  *  \ 

b2  =  T— - 5  ^  (Oj^xj  exp  [-  ico^Zj  , 

C  ^2 

C2  =  -^(fsina)^x-cos6)^x)exp(-f022)  • 

^2 

The  quantities  (bi,  Ci)  and  (&2/  ^2)  be  expressed  in  terms  of  incoming 
and  outgoing  waves  as  in 


_1 

2  co^(o^2(oli 


[{(0x2  -  (Oxl)  +  («x2  +  (Oxi)  exp  (-  ico^z)  , 


- A{f^l(0ja  -  f^2(0xl]  +  {f^l(0x2  +  f(^2(0xl]  exp  (-  ico^z)  , 

iK^CO^COxi 


&2  =  -T 


K2(gj/c' 


(O^col2 


e  ““*2*  exp  (- lODjZ)  , 


62  =  — —  e  ’"^2^  exp  (-  ico^z)  . 
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Because  of  the  lack  of  dependence  of  {b,  e)  on  y,  the  equations  for  the  elec¬ 
tric  and  magnetic  fields  given  in  equation  (29)  simplify  to 


D  ^ 

^x-CO,^  ' 

r  ^ 

^x-CO,^  ' 

u  KO)de 

„  ab 

B,  =  ik^b  , 

E,  =  ik\  , 

(55) 


where  {k,  k^)  stand  for  (Ki,  k^)  or  (k2,  k^),  depending  on  the  medium  in 
which  the  fields  are  evaluated.  Equations  (53)  and  (54)  are  substituted  into 
equation  (55)  to  produce  values  of  the  B  and  E  fields. 

The  parameter  T  defines  the  polarization  modes.  If  T  is  large,  e,  by  com¬ 
parison,  is  small,  and  By,  E^,  and  are,  effectively,  zero.  This  condition 
defines  the  mode  in  which  the  electric  field  is  perpendicular  to  the  plane  of 
incidence  and  refraction.  Apart  from  factors  that  are  common  to  both 
regions,  the  fields  for  this  mode  in  region  1  are 


=  co^i] +  ^(o, +  0),^) , 

Bz  =  -  ^  (^x2  -  -  \  CO,,  {co^  +  co„)  , 

Ey  =  I  -  co„)  - 1 0){o)^  +  (o„)  , 


and,  in  region  2,  are 


Ey  =  -(oo},2e~^^^^’‘  . 


(57) 


Likewise,  if  T  is  zero,  b  is  zero  in  both  media  and  B„  B,,  and  Ey  are  zero. 
This  value  of  T  defines  the  polarization  mode  in  which  the  magnetic  field 
is  perpendicular  to  the  plane  of  incidence  and  refraction.  The  fields  for  this 
mode  in  region  1  are 


®y  =  ^  7  ('"1  ^X2 -  ^^2  co„)  -^^{k,co,2  +  K2  (»„)  , 


2  c 
1  (0^ 


2  c 

1  co^, 


=  2  —  ('^1  ^x2 -  ^^2  CO,,}  +  ^2  0),,}  e 


-IO),,X 


1  CO„C 


^-^(>^1^x2-  ^2  COxl)  ^  (k-1  ®x2  +  X’2  «xl)  ^  " 


(58) 


-t6)x,X 
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and,  in  region  2,  are 


=  (59) 

The  Fresnel  equations  are  usually  expressed  as  ratios  of  reflected  and 
refracted  fields  to  the  incident  fields,  and  the  fields  are  expressed  in 
terms  of  the  angle  between  the  direction  of  propagation  and  the  surface 
normal.  The  geometry  is  illustrated  in  figure  4,  which  shows  the 
wavefronts,  the  direction  of  propagation,  and  the  coordinate  system  in 
which  the  x-coordinate  is  normal  to  the  surface.  The  spatial  frequencies 
and  cOj  may  be  expressed  in  terms  of  the  angle  0as 


(0^ 


(60) 


Substitution  of  these  relationships  into  equations  (56)  through  (59)  leads  to 
the  more  common  form  of  the  Fresnel  equations  [15]. 


Figure  4.  Definition  of 
angle  6  that  direction 
of  propagation  makes 
with  surface  normal. 
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6.  Boundary  Conditions  at  a  Conducting  Surface 

The  first  group  of  solutions  for  boundary  conditions  at  a  dielectric  inter¬ 
face  (eq  (17))  can  be  used  to  describe  boundary  conditions  at  the  surface  of 
a  conductor;  these  potentials  satisfy  the  boundary  conditions  of  equation 
(11).  Substitution  of  equation  (17)  into  equation  (12)  lead  to  linear  equa¬ 
tions  for  the  constants 


{<l>-(od)^-Q)a^=-£-  , 

dx  dy  Eo 

I  jsdu  do 

-ik^av  =  -njy  , 

where  the  subscripts  on  (^,  a^,  d,  a)  that  denote  the  region  and  the  set  of 
solutions  have  been  dropped.  The  Lorentz  condition  is  invoked,  as  before, 
to  eliminate  the  constant  d,  leading  to  three  equations  for  the  three  remain¬ 
ing  constants.  In  matrix  form,  these  equations  are 


QXO^ 

k^ 

QXj 

<l> 

p/£o 

du 

0)1(0/C^) 

dx 

0)^ 

Piz 

/ 

(62) 

0 

0 

•7  2  V 

3u/a«| 

a 

Piy 

where,  as  before,  q 
matrix  equation  is 


do/dy 

du/dx 


A  set  of  linearly  independent  solutions  to  the 


(p 

co/c 

a)z 

«z 

(Ojc 

f 

fo/c^ 

/ 

(olq/o) 

a 

0 

0 

-[coco^yickf 

The  first  solution  has  an  eigenvalue  of  zero  and  defines  gauge  transforma¬ 
tions.  The  third  solution  differs  from  that  of  equation  (21).  It  will  be  shown 
that  these  solutions  define  TM  and  TE  propagation  modes,  where  the 
direction  of  propagation  is  in  the  z-direction.  The  modes  resulting  from 
these  solutions  will  be  the  same  as  those  modes  usually  employed  to  ana¬ 
lyze  conducting  waveguides;  however,  any  pair  of  linearly  independent 
solutions  that  lead  to  nonzero  charges  and  currents,  and  that  is  also  inde¬ 
pendent  of  the  gauge  solution,  defines  physical  field  configurations. 
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The  electric  and  magnetic  fields  are  computed  as  before  by  substituting 
equation  (63)  into  equation  (28).  The  fields  are 


and 


0)  ^zdv 

' 


_  or  0)z  do 


Although  the  vector  potentials  were  computed  at  a  single  point  of  the 
boundary,  the  only  quantity  that  characterizes  individual  points  (the  ratio 
q)  does  not  appear  in  the  expressions  for  the  fields.  Unlike  the  dielectric 
interface,  computations  at  any  point  on  the  conducting  surface  lead  to  the 
same  equations  for  the  fields,  even  if  q  is  not  constant.  Moreover,  as  a  con¬ 
sequence  of  the  rotationally  invariant  forms  in  equation  (17),  the  forms  of 
equations  (66)  and  (67)  are  also  invariant  to  rotations  in  the  {x,  y)  plane. 
These  are  the  equations  for  the  fields  in  a  conducting  waveguide  [16],  rela¬ 
tive  to  any  coordinate  system  whose  z-axis  is  in  the  direction  of  propaga¬ 
tion.  With  different  variables,  they  also  describe  other  two-dimensional 
problems. 


24 


7.  Discussion 


The  three  orthogonal  magnetic  vector  potentials  are  represented  as  [9] 


A^='Vy/i 
^,  =  Vx(iv^) 

VxVx(sv'^,]  , 

where  s  is  an  arbitrary  constant  vector  that  has  been  called  the  pilot  vector 
[17]  and  y/;-  is  a  scalar  solution  to  the  wave  equation.  The  electric  potential 
<I>  can  be  constructed  from  the  other  components  of  A^,  A^^,  or  A^  by  using 
the  Lorentz  condition  (eq  (7)).  The  vector  A^  and  the  potential  O  that  is 
associated  with  it  define  gauge  transformations;  the  electric  and  magnetic 
fields  associated  with  this  set  of  solutions  are  zero.  For  each  Fourier  com¬ 
ponent  (0^,  the  subscript  in  equation  (66)  can  be  1  or  2,  because  there  are 
two  independent  solutions  to  the  scalar  wave  equation.  Excluding  the 
gauge  solutions,  four  sets  of  solutions  to  Maxwell's  equations  are  repre¬ 
sented  in  equation  (66).  This  representation  is  consistent  with  the  four 
solutions  in  equations  (30)  and  (31).  Different  sets  of  scalar  potentials  arise 
through  choices  of  the  pilot  vector  i ,  the  independent  scalar  solutions 
and  the  choice  of  gauge. 

The  Debye  potentials  can  be  related  to  equation  (66)  by  choosing  the  pilot 
vector  s  in  the  x-direction,  normal  to  the  interface.  The  solutions  to 
Maxwell's  equations  in  this  case  are 


<I> 

mxj/i 

0 

0 

dy/i/dz 

-d¥i/^ 

d^y/j/dxdz 

Ay 

d¥i/^ 

/ 

dy/fldz 

/ 

a  y/i/dxdy 

A. 

dy/j/dx 

0 

-(aV,/5f/^  +  5V,/3fe^) 

where  y  and  z  are  directions  that  are  locally  tangent  to  the  interface.  Substi¬ 
tution  of  the  second  and  third  solutions  of  equation  (67)  into  the  equations 
for  the  electric  and  magnetic  fields  (eq  (29))  leads  to  sets  of  solutions  with 
electric  or  magnetic  field  components  of  zero  in  the  x-direction.  This  prop¬ 
erty  was  used  by  Debye  to  define  these  potentials  [5].  By  using  the  conti¬ 
nuity  of  tangential  components  of  E  and  H  at  a  dielectric  interface,  Debye 
derived  boundary  conditions  that  the  potentials  obey  at  the  surface  of  a 
sphere  of  radius  a  as 


a^rv^ii  (fcifl)  _a^rvA2i  (^2«) 
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where  the  first  subscript  refers  to  a  dielectric  region  and  the  second  sub¬ 
script  refers  to  the  scalar  solution  to  the  wave  equation.  Debye  computed 
scattering  intensities  by  constructing  solutions  that  obeyed  the  radiation 
condition  in  the  outer  region  that  were  finite  in  the  inner  region,  and  that 
obeyed  the  botmdary  conditions  of  equation  (68).  If  the  vector  s  lies  in  the 
2-direction  in  the  plane  of  the  interface,  a  set  of  solutions  is 
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If  one  of  the  scalar  potentials  is  taken  to  be  u  (x,  y)  and  the  other  v  {x,  y), 
these  solutions  can  be  combined  to  give  the  potential  forms  of  equations 
(17)  and  (22).  For  equation  (17),  the  gauge  solution  is  combined  with  the 
last  solution  (both  with  scalar  function  u  {x,  y),  and  the  combination  is 
added  to  the  second  solution,  with  scalar  function  v  {x,  y).  The  forms  of 
equation  (22)  are  found  in  a  similar  way. 

The  coupled  azimuthal  potentials  developed  by  Morgan,  et  al  [10]  can  also 
be  described  by  equation  (66).  They  were  developed  using  cylindrical 
coordinates,  although  the  development  can  be  extended  to  other  coordi¬ 
nate  systems.  In  the  notation  and  unit  system  used  here,  fields  resulting 
from  the  coupled  azimuthal  potentials  are 
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and  the  azimuthal  dependence  is  of  the  form  exp(im^),  where  m  is  an  inte¬ 
ger.  The  permittivity  e  and  the  permeability  n  are  assumed  to  be  functions 
of  the  coordinates  (r,  z).  If  they  were  constant,  and  if  the  quantity 
were  not  in  the  definition,  equation  (70),  without  this  factor,  could  be 
derived  from  potentials  of  the  form  of  equation  (66),  with  the  pilot  vector 
s  in  the  (j)  direction.  In  this  case,  the  potentials  would  not  be  coupled,  but 
would  satisfy  the  wave  equation  individually. 

Bouwkamp  and  Casimer  [18]  derived  a  relationship  between  the  Debye 
potentials  and  current  distributions.  They  showed  that,  if  the  fields  of  the 
eigenfunctions  of  a  dielectric  sphere  were  generated  by  currents,  the  fields 
resulting  from  an  electric  multipole  distribution  of  currents  would  be 
derivable  from  one  of  the  Debye  potentials,  and  those  fields  resulting  from 


magnetic  multipole  distributions  would  be  derivable  from  the  other.  The 
scalar  potentials  u  {x,  y)  and  v  (x,  y)  can  be  related  to  surface  currents.  Com¬ 
ponents  of  equation  (21)  satisfy  the  boundary  conditions  of  equation  (11) 
and  can  be  associated  with  currents.  These  components  can  be  substituted 
into  equation  (61).  The  first  component  is  the  gauge  solution  and  leads  to 
charges  and  currents  of  zero.  The  second  component  leads  to  the  equations 
for  surface  currents 


and  the  third  leads  to 
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8.  Conclusions 


The  preceding  sections  describe  an  approach  for  computing  electric  and 
magnetic  fields  using  scalar  solutions  to  the  wave  equation  that  are  zero  at 
an  interface  and  have  normal  derivatives  of  zero  at  the  interface.  The  use  of 
these  two  functions  simplifies  their  matching  at  the  interface.  The  scalar 
functions  are  computed  from  the  locations  of  the  interface  and  the  con¬ 
stants  (£,  ju)  that  characterize  the  regions.  Boundary  conditions  for  the  elec¬ 
tric  and  magnetic  fields  have  been  satisfied  by  forming  linear  combina¬ 
tions  of  pairs  of  scalar  functions  that  solve  the  wave  equation. 

The  computations  presented  here  are  two-dimensional.  At  a  dielectric 
interface,  it  is  assumed  that  the  ratio  q  of  the  derivatives  of  the  scalar  func¬ 
tions  is  constant  along  the  boundary  of  the  two-dimensional  region,  and 
that  the  normal  derivatives  of  u  are  proportional  on  each  side  of  the  inter¬ 
face.  These  conditions  imply  that  values  of  v  are  also  proportional  on  each 
side  of  the  interface.  The  coordinate  system  must  be  chosen  to  ensure  that 
these  conditions  hold.  These  conditions  replace  the  matching  formulas  of 
equation  (68). 

Many  analyses  of  boimdary  conditions  at  a  conducting  surface  have  been 
carried  out  in  terms  of  u  (x,  y)  and  v  {x,  y),  although  these  functions  are  not 
regarded  as  scalar  potentials,  but  as  representations  of  tangential  electric 
and  magnetic  fields,  since,  in  this  case,  the  tangential  fields  obey  exactly 
the  same  boundary  conditions.  It  has  been  shown  here  that  these  scalar 
potentials  can  also  be  applied  to  dielectric  interfaces,  and  that  they  can 
assist  in  simplifying  computations  in  this  case. 
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